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ABSTRACT 

We have analyzed the double-lined eclipsing binary system ASAS J180057-2333.8 from the 
All Sky Automated Survey (ASAS) catalogue . We measure absolute physical and orbital pa¬ 
rameters for this system based on archival V-band and /-band ASAS photometry, as well as 
on high-resolution spectroscopic data obtained with ESO 3.6m/HARPS and CORALIE spec¬ 
trographs. The physical and orbital parameters of the system were derived with an accuracy 
of about 0.5 - 3%. The system is a very rare configuration of two bright well-detached giants 
of spectral types K1 and K4 and luminosity class II. The radii of the stars are Ii\ = 52.12 ± 
1.38 and R 2 = 67.63 ± 1.40 R 0 and their masses are Mi = 4.914 ± 0.021 and M 2 = 4.875± 
0.021 M 0 . The exquisite accuracy of 0.5% obtained for the masses of the components is 
one of the best mass determinations for giants. We derived a precise distance to the system of 
2.14 ± 0.06 kpc (stat.) ± 0.05 (syst.) which places the star in the Sagittarius-Carina arm. The 
Galactic rotational velocity of the star is 0 S = 258 ± 26 km s -1 assuming ©o = 238 km s -1 . 
A comparison with PARSEC isochrones places the system at the early phase of core helium 
burning with an age of slightly larger than 100 million years. The effect of overshooting on 
stellar evolutionary tracks was explored using the MESA star code. 

Key words: binaries: eclipsing - binaries: spectroscopic - stars: fundamental parameters - 
stars: individual: ASAS J180057-2333.8. 


1 INTRODUCTION 

The analysis of binary star systems is a very important part of astro¬ 
physics. Calculations of their orbits allow us to directly determine 
the masses of their components, which also gives us a chance to es¬ 
timate other physical parameters. Moreover, analysis of eclipsing 
binary systems can provide the absolute values for the numerous 
physical parameters of stars, which are essential for testing stellar 
structure and evolutionary models. In particular, the SB2 eclips¬ 
ing binary systems allow us to directly deter mine those parameters 
with a requisite accuracy dTorres et alj2010l) . From both photomet¬ 
ric and spectroscopic data, we can measure very accurate distance- 
independent stellar parameters such as stellar masses, radii, lumi¬ 
nosities and effective temperatures. 


* E-mail: ksenia@astrouw.edu.pl 


In this paper we present the first determination of the physical 
and orbital parameters of the double-lined detached eclipsing bi¬ 
nary system from the All Sky Automated Survey (ASAS ) identified 
as AS AS J180057-2333.8 (hereafter ASAS1800) by IPoimanskil 
d2002l) . The star is also cataloged as TYC 6842-1399-1 and 2MASS 
J18005707-2333420, and classified as a detached binary system in 
the AC VS catalog ue. Its V magnitude at maximum brightness is 
10.19 dPoimanskil 12002 3 and the amplitude of photometric vari¬ 
ations in the V-band is 0.47 mag. It has a circular orbit with a 
period of 269 days. The system contains two evolved giants. It 
is located in the Galactic disc (b = —0°2) and has not been 
spectrosco pically analyzed before . Such systems are very rarely 
found (eg. iHelminiak et al.1 l2015h in our Galaxy. Although it is 
relatively young, the binary system is composed of well detached 
bright giant stars. Late type eclipsing binary systems are one of 
the best candidates for distance determinations dPietrzvnski et al.l 
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120131 : [Thompson et al]|200~i1 ). The precision of our distance deter¬ 
mination to this eclipsing binary (total error of ~ 4%) rivals those 
obtained from interferometric parallaxes of Galactic masers (e.g. 
Xu et al] |20li their Table 1) at comparable distances. In this pa¬ 
per we focus on a precise determination of physical parameters of 
the system, its distance and space kinematic properties and a dis¬ 
cussion of evolutionary status. Absolute dimensions of both com¬ 
ponents are used in discussion of ASAS1800’s evolutionary status. 
We begin with a presentation of the data collection and analysis fol¬ 
lowed by a description of our results. In the last section we present 
our conclusions. 


2 OBSERVATIONS 
2.1 Photometry 

For our analysis of ASAS1800 we used the arc hival V-band and 
/-band photometry from the ACVS |Poimariskill200Cf) . A total of 
887 and 266 measurements were obtained in the V-band and the 
/-band, respectively, and the data coverage for the light curve for 
this system is complete in both filters. The primary eclipse is to¬ 
tal. The time span of the observations for the V-band is 3189 days 
(JD 2451949 to JD 2455138) and 1675 days for the I -band (ID 
2452282 to ID 2453957). The magnitude in the K-band was taken 
from 2MASS catalogue and is K = 5.917 mag dCutri et aill2003l ). 
The observation was made at an orbital phase of <f> = 0.125, well 
separated from either eclipse. 


2.2 Spectroscopy 

The high-resolution spectra were collected with the ESO 3.6 tele¬ 
scope at La Silla Observatory, Chile equipped with the HARPS 
spectrograph, as well as with the Euler 1.2m telescope at La Silla, 
Chile, equipped with the CORALIE spectrograph. The resolution 
of the CORALIE spectrograph is ~50,000. The HARPS spectro¬ 
graph was used in the EGGS mode at a resolution of ~80,000. For 
our analysis we used 14 spectra in total, 12 of which were taken 
with the HARPS spectrograph and 2 spectra with CORALIE. 


3 ANALYSIS AND RESULTS 


In order to derive absolute physical and orbital parameters for the 
system, we used the Wilson-Devinney code (WD) , version 2007 
1 van Hamme & WilsorJ 120071 : IWilson & Devinnevl 1 197 ll : Iwilsonl 
ll979L 1990lf . equipped with the automated differential correction 
(DC) optimizing subroutine and Monte Carlo simulation package. 
The WD code allows us to simultaneously solve multiband light 
curves and radial velocities which is recommended as the best way 
to obtain a consistent model of a binary sy stem (Wilson 2007) . We 
also used RaVeSpAn software written by IPilecki et alj l2012t) for 
measuring radial velocities, as well as for spectrum disentangling, 
used in further analysis. 


3.1 Radial Velocities 

The RaVeSpAn code uses the Broadening Function formalism 
dRucinskilfpwil 19991) to measure radial velocities of the compo¬ 
nents of the binary. Tem plates w ere se lected from the synthetic 
library of LTE spectra of ICoelho et al.l d2005i) . We calculated the 
components’ radial velocities over the wavelength range 4360 to 


Table 1. Radial velocity measurements. Typical uncertainty is 40 m s 1 . 


HJD 

-2450000 

Vi 

(km s -1 ) 

v 2 

(km s -1 ) 

Instrument 

5448.48918 

-36.806 

1.770 

HARPS 

5449.55107 

-37.554 

2.520 

HARPS 

5467.51220 

-47.416 

12.664 

HARPS 

5470.49095 

-48.694 

13.910 

HARPS 

5478.53899 

-51.214 

16.369 

HARPS 

5479.51462 

-51.498 

16.601 

HARPS 

5499.50211 

-51.982 

17.304 

CORALIE 

5500.49960 

-51.831 

17.005 

CORALIE 

5502.51109 

-51.492 

16.653 

HARPS 

6214.50228 

-3.458 

-31.415 

HARPS 

6448.68949 

15.524 

-50.835 

HARPS 

6449.83022 

15.196 

-50.472 

HARPS 

6450.84722 

14.856 

-50.261 

HARPS 

6553.63729 

-50.404 

15.572 

HARPS 


6800 A, excluding atmospheric and strong hydrogen lines. The re¬ 
sulting radial velocity curve is presented in Fig. |T| and the mea¬ 
sured radial velocities are presented in the Tab. Q] The components 
differ in systemic velocity, and we applied a correction of v sys = 
-175 m s _1 s to the radial velocities of the secondary component 
in order to obtain an accurate radial velocity solution. Such a shift 
can be caused either by a differential gravitational redshift between 
the stars or due to large-scale convective motions (eg. iTorres et~ahl 
120091) . 


We also determined the rotational velocities of both compo¬ 
nents in order to compare them with the expected synchronous 
velocities. We determined rotational velocities using the RaVeS¬ 
pAn code, fitting to rotationally broadened profiles. The measured 
broadenings are vm 1 = 11.02± 0.16 km s^ 1 and vm 2 = 14.84 ± 
0.28 km s -1 . To determine v sin* we also had to take into account 
the macroturbulence and instrumental profile contribution to our 
earlier measurements. In order to estim ate those values we used the 
rela tion presented in Massarotti et al.l d2008t see their Equation 1) 
and lTakeda et al .1 d2008tl (» m i = 0.42/^). The measured velocity 
is linked with rotational velocity through the relation: 


Vm — Vrot + (0.42£ht) 2 + Dip 


( 1 ) 


where £rt is the radial-tangential macroturbulence, and Vi p is the 
instrumental profile broadening. We estimated the macroturbulence 
to be Vmti = 3.20 km s -1 and Vmt2 = 2.61 km s -1 . The instrumen¬ 
tal profile was estimated to be Vi p = 2.25 km s -1 , assuming the 
resolution of the HARPS spectrograph i n EGGS mode to be R = 80 
000 and using a relation given in lTakeda et al .1 J2008h . With these 
assumptions, we estimated the rotational velocities to be Vi sin i = 
10.31 ±1.16 and V 2 sin i = 14.44 ± 1.28 km s -1 . 

Assuming that the rotation axes are perpendicular to the or¬ 
bit and the rotation is synchronized with orbital motion, the ex¬ 
pected equatorial, rotational velocities are vi = 9.79 km s^ 1 and V 2 
= 12.71 km s" 1 . We calculated v using a formula: 


2ti -R 


( 2 ) 


where R is the radius of a component and P is the orbital period of 
the system. The measured rotational velocities are consistent with 
the expected synchronous velocities within 0.5a and 1.4a for the 
primary and secondary components, respectively. We conclude that 
the components are rotating synchronously. 
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Orbital phase 


Figure 1. Radial velocity curve solution to ASAS 1800 from the WD code. 
Filled black circles represent measurements of the primary, the empty cir¬ 
cles measurements of the secondary. The empty triangle stands for the mea¬ 
surement of the secondary, observed at HJD 2456214.50228, which was not 
taken into account. 


Table 2. Atmospheric parameters of the components. 


Component 

Teff IK] 

[Fe/H] 

logs 

Primary 

4535 ± 70 

-0.14± 0.1 

1.88 

Secondary 

4240 ± 70 

-0.27 ±0.1 

1.93 


3.2 Spectral disentangling and atmospheric analysis 


The spectral disentangling was done using a method outlined by 
iGonzalez & Levatg j 2006 ). We u sed the two step method described 
in detail in iGraczvk et al di2014 T to derive properly renormalized 
disentangled spectra. These spectra were used for deriving the ba¬ 
sic atmospheric parameters of effective temperature T e s, gravity 
logg, microturbulance vt, and metallicity [Fe/H] assuming the 
local thermod ynamical equilibrium and using pro gram MOOG 
i Snedenl ll973h . Details o f the method ar e given in iMarino et al .1 
(2008) and the line list in I Villanova et al] feOldt . Values of these 
parameters are presented in Tab. [2] The derived effective tempera¬ 
tures of both components were used in the WD model and to cal¬ 
culate interstellar reddening. 

Taking into account the metallicity determination uncertainty 
and additional inaccuracies connected with spectra disentangling 
and renormalization, the difference in the metallicities of both com¬ 
ponents can be neglected (Tab. [2^. Therefore, we can assume that 
the components have common metallicity within the margin of er¬ 
ror. In our analysis of the evolutionary status of ASAS 1800 we as¬ 
sumed the metallicity of the system to be equal to the metallicity of 
the primary component - [Fe/H] = — 0.14 dex (see Section [X6l >. 


3.3 Interstellar extinction 


We estimated reddening b ased on several c a librations of 
T e ff - (V — K) colour ddi Benedettol 1 19981: lAlonso et al.1 


1999); Houdashelt et alj 2000J; Ramirez & Melendez 

2005; 

Gonzalez Hernandez & Bonifacio 20091: Casssrande et al 

2010; 

Masana et al. 2006), using the values of T e n presented in Tab. [2] 


We determined E(B - V) = 0.525 ± 0.035mag. The errors on 
E(B-V) result from the accuracy of our effective temperature 
determination (see Tab. [3 and front the accuracy of the adopted 
effective temperature - colour calibrations, color 

We also determined the interstellar extinction from from the 


calibration of effective temperature - (V — I) colour. From the ef¬ 
fective temperature - {V — I) colour calibrations of lWorthev & Lee! 
d20 1 ill we estimated the intrinsic (V — I) colours of each compo¬ 
nent. These colours were then compared with the observed colours 
of the components to estimate the E(V -1) extinction. We estimated 
the reddening E(V-I) for both components and then transformed it 
to E(B-V) using: 


E(B -V) = 


E(V - I) 


(3) 


1.399 

The mean value of the reddening was E(B-V)=0.609 ± 0.042 . 

Finally, we use d extinctio n map_s_of Schlegel et ~al] d 1998h with 
the recalibration of Ischlaflv & Finkbeiner ! (2011 ) to estimate the 
reddening in the direction of ASAS 1800. The total foreground red¬ 
dening in this direction is E(B - V) = 26.597 mag. Since the red¬ 
dening to ASAS 1800 is only a fraction of this number, we have to 
assume a distribution of dust within the Milky Way and to know 
the distance to our system. The simple axisymetric model of the 
exponential disc gives a density of matter within the Galaxy: 


p(r,z) = po exp(—r/rd - \z\/zd) 


(4) 


where rj, and Zd are the disc scale len gth and height, respectively . 
We adopted the following values from Drimmel & Spergell d200lh : 
sun’s height in the Galactic disc ho = 0.015 kpc, rd = 3.2 kpc and Zd 
= 0.135 kpc. The Galactic coordinates of ASAS1800 are l = 6.37° 
and b = —0.23°. Moreover we assumed th e solar distance to the 
Milky Way centre to be Ro = 8.3 kpc dpi I lessen et al.ll2009h and 
the distance to ASAS1800 to be D = 2.16 kpc (Section 13.8b . We 
assumed a Milky Way disc truncation at D ou t er = 20 kpc. Writing 
r and z as functions of distance d f rom the Sun in the direc tion 
of ASAS 1800 we obtain r(d) = \J (7?q + d 2 — 2Rod cos l) and 
z(d) = \h 0 + dsin&|. Substituting those functions into Eq. [4]we 
obtain the relation p = p(d). We numerically integrated this re¬ 
lation along the line of sight twice: from 0 to D, corresponding 
to the reddening of the eclipsing binary, and from 0 to D outer , 
corresponding to the foreground reddening. The ratio gives E(B - 
V)asasi800 = 0.440 ± 0.057 mag. 

We adopt a final value of E(B-V) = 0.52 ± 0.07, the mean 
value from all our estimates. The error is a combination of both 
statistical and systematic error, dominated by the statistical error. 


3.4 Modeling 

The WD code is based on Roche lobe geometry and employs a so¬ 
phisticated treatment of stellar surface physics. It fits a geometric 
model of a detached eclipsing binary to a light curve in order to 
establish parameters of the system and its components. The orbital 
period and the moment of primary minimum were derived from the 
AC VS data. We measure P = 269.363 days and T 0 = 2452728.5. 
The moment of the primary minimum (To) was later adjusted dur¬ 
ing the further analysis. The average out-of-eclipse magnitudes 
were established taking into account all of the observational data 
outside of minima. We measure V = 10.319 mag and I = 8.231 mag. 
Since the primary eclipse is total we were also able to directly deter¬ 
mine the magnitudes of the components: Vs=l 1.097 and Is=8.935 
for the secondary, and Vp = 11.061 and Ip = 9.071 for the pri¬ 
mary. We refer to the primary component as the star which is being 
eclipsed in the deeper, primary minimum. 

We simultaneously fitted two light curves, in the J-band and 
V -band, as well as the radial velocity curves. The in put parameters 
for the DC subroutine were chosen as described in IGraczvk et al .1 
j2012l) . When using the Wilson-Devinney code, it is important to 
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carefully define which parameters are adjustable in order to arrive 
at the best fitted model. In our analysis we decided to adjust the or¬ 
bital semi-major axis (a), systematic radial velocity ( 7 ), the orbital 
inclination (i), the average surface temperature of the secondary 
component (T 2 ), the modified surface potential of both components 
(fix, fl 2 ), the mass ratio (q = M 2 /M 1 ), time of the primary min¬ 
imum (To), the observed orbital period ( P 0 bs ), and the relative lu¬ 
minosity of the primary component in the two bands (Lly, LI/). 

To set the effective temperature scale for each component, we 
ran the WD code with the initially assumed temperature of Ti = 
4700 K assuming a spectral type of K0 III. From the preliminary so¬ 
lutions obtained from the WD code we derived approximate surface 
gravities for the components of the binary of log g\ = 1.7 and log 
<72 = 1.4, as well as luminosity ratios in the V, I and TT-bands. The 
resulting luminosity ratios together with reddening E(B-V) (Sec¬ 
tion 13.31 were used to obtain the dereddened (V — K ) colour 
index. The 2MASS magnitudes were converted onto the Johnson 
photometr i c syst em usi ng updated tran s forma tion equations from 
ICarnenteil 1200 llFI and iBessell & Brettl d 1988h (Tab. [ 6 }. Knowing 
the approximate log <?i = 1.7, the dereddened {V — K) = 2.61 and 
assuming [Fe/H] = 0 dex, we were able to estimate preliminary ef¬ 
fective t emperatures of the com ponents based on the calibrations 
given by IWorthev & Leej d201 lh . The resulting effective tempera¬ 
ture was set as the temperature of the primary Ti = 4550 K. We 
used this value as a starting point for our analysis and then iterated 
to find the best solution for both the V -band and /-band light curves 
using the LC subroutine of the WD code. All free parameters were 
adjusted at the same time. 

The albedo and gravity brightening parameters were set to 
0.5 and 0.32 respec tively, which are appropriate values for those 
kind of stars ILucvll 1 9671) . To compute the lim b dar kening co ef¬ 
fi cient s we used the logarithmic law of iKlinglesmith & Sobieskil 
dl97C)h . Those coefficients were calculated internally by the WD 
code during each i teration of DC using tabulated data computed by 
Ivan Hammel d 19931 ). Additionally, we calculated models using the 
linear and square root limb darkening law. However, that resulted in 
a slightly worse fit to the light curves, changing the stellar parame¬ 
ters of the system by less than 0.5%. Thus, we adopted the solution 
obta ined with fixed coeffici ents of the logarithmic limb darkening 
law jPietrzvnski et al.ll2013l ). 

At the end of the fitting procedure we additionally adjusted the 
third light (J 3 ) to determine its impact on the solution. Formally, the 
solution suggested an unphysical value for I 3 , and we therefore set 
I 3 = 0 in our final solution. 

The solution, especially the luminosity ratio of the compo¬ 
nents, was used to renormalize the disentangled spectra. Subse¬ 
quently, the atmospheric analysis was performed in order to obtain 
a better estimation of the temperatures of the components and their 
metallicities (see Section [372). We derived effective temperatures 
of Ti = 4535 ± 70 K and T 2 = 4240 ± 70 K and metallicities of 
[Fe/H] 1 = -0.14 ± 0.1 dex and [Fe/H ] 2 = -0.27 ± 0.1 dex. We then 
adopted Ti as the new effective temperature of the primary com¬ 
ponent and we repeated the fitting using the DC subroutine of the 
WD code. 

The /-band light curve solution obtained with the WD code 
is presented in Figure [2] and for the V -band in Figure [3] and the 
parameters are summarized in Table [3] 


0 0.2 0.4 0.6 0.8 1 



Figure 2. The /-band light curve of ASAS1800 together withe the solution 
from the WD code. 


0 0.2 0.4 0.6 0.8 1 



Figure 3. The V -band light curve of AS AS 1800 together with the solution 
from the WD code. 


Table 3. Photometric and orbital parameters obtained with the Wilson- 
Devinney code. 


Parameter 

WD result 

Orbital inclination i (deg) 

88.67 ± 0.21 

Orbital eccentricity e 

0.0 (fixed) 

Sec. temperature T 2 (K) 

4211 ± 13 

Fractional radius r\ 

0.1387 ±0.0020 

Fractional radius 7*2 

0.1800 ±0.0012 

(n ± r 2 ) 

0.3187 ±0.0012 

k = r 2 /n 

1.2976 ±0.0104 

Observed period P Q bs (day) 

269.496 ±0.014 

(L2/L1) V 

0.9850 ± 0.0054 

(L2/L1)/ 

1.1650 ±0.0072 

(L2/Ll)j 

1.3379 

(L2/L1) k 

1.5249 

T 0 (JD-2450000) 

2728.82 ± 0.06 

Semimajor axis a (Rq) 

375.72 ± 0.37 

Systemic velocity 7 (km s — *) 

-17.625 ±0.021 

Prim, velocity semi-amplitude K i(km s — - 1 ) 

35.11 ±0.10 

Sec. velocity semi-amplitude if 2 (km s —1 ) 

35.38 ±0.10 

Mass ratio q 

0.992 ± 0.003 

RV rmsi (km s — - 1 ) 

0.056 

RV rms 2 (km s — - 1 ) 

0.042 


1 http://www.astro.caltech.edu/~jmc/2mass/v3/ 
transformations/ 
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Table 5. Error budget of the distance moduli of the ASAS 1800. 


Type of error 

(m — M) 

crA 

cr(MonteCarlo) 

a diBenedetto 

aE{B - V) 

aV 

oK 

(L 2 /L!)k 

Combined Error 


(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

Statistical 

11.655 

0.003 

0.049 

- 

0.024 1 

0.0241 

0.019 

- 

0.063 

Systematic 

11.655 

- 

- 

0.043 

- 

0.03 

0.01 

0.053 


1 - combination of statistical and systematic error 


Table 4. Physical Properties of the ASAS 1800. 


Property 

The Primary 

The Secondary 

Spectral type 

K1 II 

K4II 

V a (mag) 

11.061 

11.098 

V — I a (mag) 

1.991 

2.162 

V — K a (mag) 

4.124 

4.598 

J — K a (mag) 

1.087 

1.229 

Radius ( Rq ) 

52.12± 1.38 

67.63 ± 1.40 

Mass (Mg) 

4.914 ±0.021 

4.875 ± 0.021 

log g (cgs) 

1.696 ±0.023 

1.466 ±0.018 

T eff (K) 

4535 6 ± 80 

421 l c ± 80 

v sin i (km s —1 ) 

10.31 ± 1.16 

14.44 ± 1.28 

Luminosity (Lq ) 

1031 ±91 

1290 ± 111 

M bo l (mag) 

-2.80 

-3.05 

M v (mag) 

-2.33 

-2.32 

[Fe/H ] 6 

-0.14 ±0.1 

-0.27 ±0.1 

E(B-V) 

Distance (pc) 

0.525 ± 0.07 

2142.5 ±63.5 (stat.) 

± 53.3 (syst.) 


b—atmospheric analysis 
c—WD solution 


3.5 Absolute Dimensions 

Table 0 gives astrophysical data about the two components. The 
physical radii of the stars result from the relation: R = r ■ a, where 
r is the fractional radius listed in Table [3] The masses are derived 
from the equations: 

Mr [MO] = 1.32068 • 10~ 2 - 1 (5) 

1 + q P- [a\ 

M 2 [M©] = Mi • q (6) 

where a is the semi major axis, q is the mass ratio and P - real 
period. The observed individual magnitudes in both V and /-band 
and V — I colour relation were derived directly from the flat bottom 
minimum of the second a ry com ponent. We used bolometric correc¬ 
tions from lAlonso et al.l i l 9991) to convert R-band magnitudes into 
bolometric magnitudes. 


3.6 Evolutionary status of ASAS1800 

In this Section we compare the physical parameters of ASAS 1800 
(Tab. 0 with results of stellar evolution calculations. We assume 
that the components of the system have common metallicity, equal 
to the metallicity of the primary (—0.14), or possibly la higher 
(—0.04). As we show below, lower metallicities lead to serious 
disagreement between the models and observations. In this ini¬ 
tial study we also assume that the masses of ASAS 1800, as de¬ 
termined in Tab. 0 are exact. In Fig. 0 we plot the PARSEC 
isochrones feressan et alJl201~3 ) for the two metallicities consid¬ 
ered. The isochrone was selected to minimize the \ 2 function in¬ 
cluding luminosities, effective temperatures and radii of the two 
components. Model values were calculated at the mass points cor¬ 


responding to the masses of the ASAS 1800 components (filled cir¬ 
cle and filled square for primary/secondary in Fig. 0. 

The comparison with the PARSEC isochrones places the sys¬ 
tem at the early phase of core helium burning. The agreement with 
isochrones is not satisfactory, however. It is better for higher metal¬ 
licity, but still the primary component is under-luminous and the 
secondary component is too cool and/or under-luminous. We note, 
however, that the PARSEC isochrones are available for only one 
fixed set of overshooting parameters, which strongly affect the evo¬ 
lutionary calculations. Below, we express the extent of overshoot¬ 
ing as a fraction of the local pressure scale height, f3 x H p . Both 
the overshooting from the hydrogen-burning core during main se¬ 
quence evolution (/3h) and overshooting from the convective en¬ 
velope (/3 e nv) affect the extent a nd the luminosity of the helium 
burning loops dAlongi et aljl99lj ). In PARSEC models, in the mass 
range considered, these are fixed at /3h = 0.5 and /3 en v = 0.7 
across the border of the convective zone determined with the 
Schwarzschild criterion. In the calculations described below, the 
extent of overshooting is measured above/below the border of the 
convective region, which is a more common approach. The re¬ 
sulting overshoot parameters roughly correspond to half of those 
adopt ed in PARSEC (see discussion in section 2.6 in lBressan et al.l 

l 2012 h . 

To explore the effect of the overshooting on stellar evolution¬ 
ary tracks we used MES A Star - a publicly availa ble stellar evo¬ 
lution code (release 6208; P axto n et alj |201 lLl2013h . Details of the 
code setup will be described elsewhere (Smolec et al. in prep.), 
here we summarize the most i mportant set tings. We use OPAL 
opacities and adopt the lAsplund et al .1 d2009l) solar distribution of 
heavy ele ments. C o nvection is modelled with the mixing length 
formalism dB ohm- Vitensel[19581) with the mixing length parame¬ 
ter resulting from calibration of the solar model (a = 1.78). The 
convective boundaries are determined with the Schwarzschild cri¬ 
terion. We account for the overshooting above the border of hy¬ 
drogen burning core, above the border of helium burning core 
(/ 3 He = 0 . 01 , fixed), and below the border of the convective en¬ 
velope. We neglect rotation, element diffusion (except in solar cal¬ 
ibration), and mass loss. For each component of ASAS 1800 we fix 
the mass (Tab. 0 and compute the evolution from the pre-main se¬ 
quence until the late AGB phase. Our small model grid consists 
of two metallicity values, —0.14 and —0.04, eight values of j3 h, 
/3h G [0.1, 0.12, 0.14, 0.16, 0.18, 0.20, 0.22, 0.24], and two 
values of jS env , /3 en v G [0., 0.35]. Along each pair of tracks (for 
the two components) we determined the models that minimize \ 2 
function including effective temperatures, luminosities and radii of 
the two components, at the same age. We first assume that the two 
components experienced the same extent of mixing at the edge 
of hydrogen burning core during their evolution, and hence have 
the same value of /3h, and then allow for a difference in /3h • In 
Figs. 0 and 0 we show our best solutions for the described two 
assumptions. In the four panels of these Figures we show the mod¬ 
els with (top) and without (bottom) overshooting from the convec¬ 
tive envelope (/3 e nv = 0.35 or /3 env = 0.) and adopting lower 
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Table 6. Out of eclipse magnitudes of AS AS 1800. 


V (mag) 

I (mag) 

J (mag) 

K (mag) 

ASAS1800 10.319 ±0.024 

Reference this work (ASAS) 

8.231 ±0.015 
this work (ASAS) 

7.178 ±0.030 i 

Cutri et al. (2003) (2MASS) 

5.942 ± 0.030 1 

Cutri et al. (.2003) (2MASS) 


1 - transformed to Johnson photometric system using equations from Bessell & Brett (T98^) and^argenter (20(H) 


(left) and higher (right) metallicity values ([Fe/H] = —0.14 or 
[Fe/H] = -0.04). 

We first analyze the models assuming the same values of 
/3h for the primary and secondary. Fig. [5] For higher metallic¬ 
ity (right panels) the helium burning loops become less luminous 
and the tracks shift toward lower effective temperatures. Conse¬ 
quently, the higher metallicity ([Fe/H] = —0.04), together with 
the smaller extent of overshooting from the hydrogen-burning core, 
mitigates the problems of an under-luminous primary and a too 
cool/under-luminous secondary, noted in the analysis of the PAR¬ 
SEC isochrones. The inclusion of the envelope overshoot in the 
models has two apparent effects on the tracks. It increases the 
vertical extent of the loops and decreases their overall luminos¬ 
ity. Hence, in the models including envelope overshooting, a larger 
extent of overshooting at the hydrogen-burning core is possible. 
We note that the matter that overshoots the convective envelope 
boundary faces a stabilizing stratification gradient. Whether the 
significant mixing is possible in such case is a subject of debate 
feressan et al.l2012l : |Pietrinferni et al.l2004)) . Modelling of evolved 
binary systems offers the best opportunity to test mixing scenarios 
at the bottom of the convective envelope. 

The best models displayed in Fig. [5] are those with over¬ 
shooting from the convective envelope (f3 e nv = 0.35) and have 
/3a = 0.18 (Fig. 0 top). The solutions assuming /3a = 0.16 are 
only slightly worse. When envelope overshoot is neglected, lower 
values of /3a (~ 0.14—0.16) are necessary. Clearly, the best models 
have higher metallicity. The inferred system’s age, given in Fig. [5] 
is very similar for all models, log(age) is within a narrow range 
from 8.007 to 8.021. As expected, the a ge is slig htly larger for 
higher metallicity models (see e.g. lSalaris & Cassini! 1 1200a) ). Also, 
larger the extent of overshooting from the hydrogen burning core, 
longer the main sequence evolution. 

We note that for the best models in Fig. [5] we nearly match 
the luminosity of the secondary, but that the primary component is 
still before the observed position. Hence, we can get a much better 
agreement between the models and observations assuming different 
values of overshooting from the hydrogen burning core during the 
main sequence evolution. A lower value of /3a for the secondary 
slows down its evolution (extends the main sequence phase) and 
allows a much better match of the system with observations at the 
helium burning phase - Fig. [6] An overshooting parameter higher 
by 0.02 — 0.04 for the secondary allows very good fits. The best are 
obtained for higher metallicity models (as in Fig. [5j. This time the 
models that neglect the overshooting from the convective envelope 
seem slightly better. The inferred ages are very similar to those re¬ 
ported in Fig. [5] Our model that matches observations best (Fig. [6] 
bottom right), assumes /3a = 0.16 for the primary and /3a = 0.18 
for the secondary and places the system at the early helium burning 
phase. 

The different extent of overshooting adopted for the compo¬ 
nents of the eclipsing binary system, with nearly equal masses and 
metallicities of the components, might appear unjustified. We note, 
however, that an overshooting parameter expresses our ignorance 
about all kinds of mixing processes that may occur at the edge of 
the convective core. In particular, rotation leads to additional mix¬ 


ing, the extent and efficiency of which depend on the rotation rate. 
We see no reason to assume that the initial rotation rate was the 
same for two stars. 

We conclude that ASAS1800 is at early phases of core helium 
burning. The age of the system is slightly larger than 100 million 
years. The models favour a metallicity that is close to solar in value. 


3.7 Tidal evolution of the system 

The observations indicate that the orbit of ASAS1800 is circular 
and rotation of both components is fully synchronized with the or¬ 
bital period, which means that the memory of the initial values of 
these parameters is at present entirely forgotten. To check, if this 
fact agrees with the predictions of the tidal evolution theory of bi¬ 
nary stars, we assume that the binary evolves as an isolated system 
with conserved total mass and angular momentum. 

Thejidal circul arization and syn chronization time scales are 
dZahnll 198^ ,1 Meibom. Mathieu & Stas sun] l2006l) 


r- tf 

(*)'■ 

(7) 

^ 21Acir C g(l + q) 


( a \ 6 

(8) 

^ 6 \ sync q 2 M R 2 

\r) ’ 


where M, R and / are respectively mass, radius and moment 
of inertia of the tidally distorted component, q is the mass ratio 
(with M in the denominator), f/ is the viscous dissipation time, 
Acirc and A sy nc are constants that depend on the internal structure 
of a star. Both time scales depend on a high power of the ratio of 
star separation to stellar radius. This ratio was of the order of 10 2 
when the binary was on the main sequence (MS), resulting in both 
time scales much longer than the MS lifetime of the components. 
Thus, the mutual tidal interaction at that evolutionary stage can be 
neglected. Substantially stronger interaction is expected when both 
components moved to the red giant region. Presently, they are both 
past the red giant tip, burning helium in their cores. We note that 
the ratio of both time scales for our binary is r C i r c /rs V nc ~ 25 — 30 
for q « 1 , A c i rc ~ A sync and I « 0.15 MR 2 dRutten & Pvlvseil 
Il988l) . That means that by the time the orbit becomes circularized, 
the components already rotate synchronously. 

To follow the eccentricity change, detailed calculations of the 
circularization rate for an evolving binary are needed. Such cal¬ 
culations have been performed by several authors for different 
kinds of systems and upper limits for periods of fully circular- 
ized binaries were obtain ed. We use the data from the paper by 
IVerbunt & Phinnevl l l 19951) who calculated the limiting period val¬ 
ues for binaries composed of giants. For giant masses correspond¬ 
ing to ASAS1800 the limiting period is equal to 616 d (see their 
Table 1). Because the period of ASAS1800, equal to 269 d, is sig¬ 
nificantly shorter than that value, we can conclude that the zero 
eccentricity of its orbit is to be expected. 

This conclusion can additionally be verified by a direct esti¬ 
mate of the absolute value of r c j rc . We use to t his purpose an ap¬ 
proximation given by IVerbunt & Phinnevl l l 19951) 
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Figure 4. PARSEC isochrones for two metallicities, [Fe/H] = —0.14 (left panel) and [Fe/H] = —0.04 (right panel). Location of primary and secondary 
components is marked with circles and squares, respectively. Filled symbols refer to the fitted isochrones and the empty ones to our measurements. 



3.80 3.75 3.70 3.65 3.60 3.80 3.75 3.70 3.65 3.60 

log F eff log r eff 


Figure 5. MESA tracks computed for the primary (M = 4.914 Mq, red, solid line) and secondary (M = 4.875 Mq, blue, dashed line) components of 
ASAS 1800, assuming the same values of overshooting from the hydrogen burning core for the primary and the secondary. Models that match the observational 
constraints best (at the same age) are marked with filled circles/squares for the primary/secondary. Models in the top two panels include convective envelope 
overshoot (neglected in the bottom panels). Metallicity is equal to [Fe/H] = —0.14 (left panels) or [Fe/H] = —0.04 (right panels). 
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3.80 3.75 3.70 3.65 3.60 3.80 3.75 3.70 3.65 3.60 

log T efs log r cff 


Figure 6. The same as Fig.[5]but for models assuming different values of overshooting from the hydrogen burning core for the primary and secondary. 


“ 3 - 4/ (3io) 4,V “” Af "‘(?) s yr ' , - <9) 

We assumed q ~ 1. If we additionally assume that the char¬ 
acteristic value of R/a « 0.2 over the giant phase, the convection 
envelope mass M env ~ M, which is a good appro ximat ion for the 
first ascend giants, T e ~ ( T\ +T2)/2 and / « 1 dZahnll 19891 ). we 
obtain r C i r c ~ 10 5 years if tides on both components are taken into 
account. This is 1-1.5 order of magnitude shorter than the lifetime 
of each component of ASAS1800 in the red giant phase so the or¬ 
bit was efficiently circularized soon after the stars reached the red 
giant branch. Because the synchronization time scale is still much 
shorter, as is shown above, the rotation of both components was 
synchronized even faster. 

3.8 Distance to the system 

To derive th e distance w e followed prescriptions given in 
iGraczvk et all d2012l [ 20 T 4 I 1 . We used V'-ban d surface brig htness 
(SF) - (V — K) color calibration measured bv ldi Benedetto! d2005t) 
for Galactic late type giant stars. The angular diameter of a star can 
be estimated using the formula: 

<(>[mas] = 10 °- 2(S - mo) ( 10 ) 

where S is the surface brightness in a given band and mo is the 
dereddened magnitude of a star in this band. We can then directly 
derive the distance to the star by scaling the angular diameter: 


d[pc] = 9.2984 • -^1 ( 11 ) 

<p[mas] 

The resulting distance to ASAS1800 is d = 2142.5 ± 63.5 
(stat.) ± 53.3 (syst.) pc (Tab. [4}. The main contribution to the sta¬ 
tistical uncertainty are random errors connected with light curve 
modelling by the WD code (connected with a relatively large dis¬ 
persion of AS AS light curves) and infrared photometry errors. Thus 
there is a significant room for improvement on the derived distance 
once high accuracy photometry will be available.. The main contri¬ 
bution to the systematic error comes from the SF calibration itself. 
The total error budget is presented in Tab[5] 


3.9 Space position and velocity 

The proper motion of the star is fi a cosS = +0.9 ± 1.75 
mas yr^ 1 and fi$ = 0.67 ± 1.71 mas yr -1 and was de- 
rived as weighted mean from three c atalogues PPMXL ca talog 
dRoeser et alj|201oll. UCAC4 catalog dZacharias et al.ll2013l) and 
SPM4.0 catalog Iciirard et al.ll201 ll) . This proper motion in Galac¬ 
tic coordinates is (/x; cos b, fib) = (T.0 +2.3, — 0.5+0.7) mas yr -1 
using the prescription given bv IPoleskH d2013l) . The calculated dis¬ 
tance corresponds to a transverse velocity in Galactic coordinates 
of (10 ± 24, —5 ± 7) km s _ . To calculat e Galactic space veloc¬ 
ity c omponents we used equations given in llohnson & Soderblomi 
d 19871) and we obtained ( u,v,w ) = (—19 + 3, 8 + 24, —4 + 7) km 
s -1 . This velocity is not corrected for solar motion with respect to 
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the Local Standard of Rest (LSR). Taking into account the peculiar 
solar mo tion (U w ,V©W©) = (11.1 ± 1.0, 12.2 ± 2.0, 7.3 ± 0.5) 
km s” 1 dSchonrich et al J|20 1 (J) and the circular speed of LSR in 
the Galaxy V c = 238 ± 9 km s -1 from lSchonrichl d2012h we ob¬ 
tain Galacto-centric velocity components of ASAS 1800 in the po¬ 
sition of the sun (Ui, Vi, Wi) = (—8 ± 4, 258 ± 26, 3 ± 7) km 
s -1 , where the errors are dominated by proper motion uncertain¬ 
ties. The Galactic space position of the star with respect to the sun 
is (A', Y. Z) = (2.12, 0.26, —0.05) kpc. The Galacto-centric dis¬ 
tance of ASAS 1800 is 6.16 ± 0.40 kpc (assuming a distance to 
the G alactic center Ro = 8.28 ± 0.38 kpc from Gillessen et al .1 
[ 2009 !) and the Galacto-centric longitude is /3 = 2°5 ± 0°3, pl acing 
the star in the Saggitarius-Carina arm (e.g. ISakai et ahll2012l . their 
Figure 3). 


4 SUMMARY AND CONCLUSIONS 

We have obtained stellar parameters for the eclipsing binary ASAS 
Jl80057-2333.8. We measure a distance to the system of 2.14 ± 
0.06 (stat.) ± 0.05 (syst.) kpc. The accuracy of the distance deter¬ 
mination is 4 %, is slightly less accurate than distances obtained 
with the same method to LMC/SMC binaries. This is due to much 
higher interstellar extinction, somewhat lower quality of the pho¬ 
tometric light curve and infrared magnitudes transformations be¬ 
tween photometric systems, leading to larger errors in the absolute 
dimensions and final distance. With better photometry and with an 
improved surface brightness - colour relation it should be possible 
to measure 1.5 - 2 % distances to such individual systems. In a re¬ 
cent s eries of papers l lPietrzvnski et al.l[2009l : iGraczvk et alj|2012l . 
l2014h we have shown that precision of 3% is already routinely at¬ 
tainable for carefully selected late type eclipsing binaries. As such 
they are a useful tool to probe the structure and the kinematics of 
the Galaxy. This technique is also an important and, moreover, in¬ 
dependent way of testing the future distance and parallax determi¬ 
nations which will be made by the GAIA mission. 

Our results also demonstrate the strength of using observations 
of well detached eclipsing binary systems in the testing of stellar 
evolution theory. Several such s ystems, wi th well dete rmined phys¬ 
ical parameters, are known (eg. lPietrzvnski et~al .20131). Evolution¬ 
ary calculations of evolved stars are sensitive to many parameters 
however, and definite conclusions require a thorough study, which 
is ongoing (Smolec et al. in prep.). 
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